Prediction of Feshbach resonances from three input parameters 
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We have developed a model of Feshbach resonances in gases of ultracold alkali metal atoms 
using the ideas of multichannel quantum defect theory. Our model requires just three parameters 
describing the interactions - the singlet and triplet scattering lengths, and the long range van der 
Waals coefficient - in addition to known atomic properties. Without using any further details of 
the interactions, our approach can accurately predict the locations of resonances. It can also be 
used to find the singlet and triplet scattering lengths from measured resonance data. We apply 
our technique to ^Li-^°K and ^°K-®^Rb scattering, obtaining good agreement with experimental 
results, and with the more computationally intensive coupled channels technique. 

PACS numbers: 33.15.Fm, 34.20.Cf, 34.50.-s 



The use of magnetic fields to control and resonantly 
enhance interactions in ultracold atomic gases allows the 
creation of highly rovibrationally excited molecules 
Alongside experimental efforts, substantial work has 
been invested in the theory of Feshbach resonances and 
molecules. The coupled channels technique [2|, (sl, [J] has 
had considerable success, but can be computationally in- 
tense. Along with the rapidly growing number of exper- 
imentally relevant collision partners |5|, [gI, [tI, [sI, [qI, lld |. 
this creates a need for accurate, computationally simple 
models of Feshbach resonances. Previous models moti- 
vated by this need have been based on, for example, an 
expansion in terms of the bound states of the singlet and 
triplet Born-Oppenheimer (BO) potentials, as with the 
very successful Asymptotic Bound State Model (ABM) 
fiot , or a fixed phase at an inner radius which is matched 
to the long range evolution of the wavefunction [Tlj . 

In this paper, we show how multichannel quantum de- 
fect theory (MQDT) can be used to construct a model 
of Feshbach resonant collisions in alkali gases that is 
both physically motivated and computationally simple. 
MQDT lH, 13, [l4] has been used to study collisions in a 
range of systems (see, e.g., [15] and references therein), 
and is built around the idea of a separation of energy 
scales. This applies to the scattering of ultracold alkali 
metal atoms because the hyperfine splitting and collision 
energies are much smaller than the potential depths. We 
will show that this enables a complete model of Feshbach 
resonances to be constructed using only the scattering 
lengths of the singlet and triplet BO potentials, and the 
common van der Waals coefficient of their long range tail. 
This represents a substantial simplification over having 
to make use of whole potentials, yet leaves us with an 
approach powerful enough to make useful predictions. In 
addition, our method is computationally simple enough 
for these three parameters to be optimized, or even found, 
from measured resonance data. 

A colliding pair of atoms can be described in terms of 



channels defined by the state of each atom, as well 
as the partial wave of the collision, ^, and its projection 
onto a quantization axis, m^. The channel energy 
is given by the energy of the two atoms at asymptoti- 
cally large separation. Typically, multiple channels are 
involved, each with its own associated potential and cou- 
plings to other channels. We use van der Waals potentials 
of the form Kdw(^) = + hH{£ + l)/{2fir^) - Ce/r^, 
where the second term represents the centrifugal barrier, 
Ce is the van der Waals coefficient, /i is the reduced mass, 
and r is the interparticle separation. This approximation 
is physically motivated by the common long range form 
of the BO potentials, which converge to the form —Cq/t^ 
beyond a certain radius r* . At r* , the energy scale of the 
potential is much greater than the channel and collision 
energies. The physics occurring in the short range region 
r < r* can then be accounted for by imposing appropri- 
ate boundary conditions at r*. In fact, the separation 
of energy scales, and accompanying length scale separa- 
tion, allow us to make the approximation r* ^ and use 
Kdw('') at all r. 

We calculate a pair of linearly independent reference 
functions, / and from each potential by solving the 
single channel radial Schrodinger equation for a chosen 
total energy, E. A channel \a) is described as open when 
E > Ea^ and closed when E < E^. An observation 
of quantum defect theory is that / and g can be ap- 
proximated to be independent of energy and angular mo- 
mentum at short range. They can be calculated numeri- 

vdw(^)' 

analytic 
approach. 

We next note that, at short range, the splitting be- 
tween the BO potentials is also far greater than the hy- 
perfine, Zeeman, and collision energies. We can therefore 
neglect these three contributions in this region. Conse- 
quently, the radial multichannel scattering wavefunction 
can be written in the form '0(r) = /(r) — K^*^^(r), where 



mentum at snort range, iney can oe caicuiatea ni 
cally , [13] or, for sonie potentials such as V^^ 
found analytically [li, 0, [111, |22[. We use the ai 
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/(r) and g{r) are vectors containing the reference func- 
tions for each channel, and bold font indicates a matrix. 
The short range K matrix, K*^*\ is to a good approxi- 
mation independent of r and E. 

Neglecting the hyperfine and Zeeman interactions at 
short range also makes K^*^ diagonal in a basis in which 
the BO potentials are diagonal. We refer to the K ma- 
trix expressed in this basis as K^^^\ For alkali atoms in 
the electronic ground state there are several ways to 
construct such a basis. A convenient choice is to first cou- 
ple together the electron spins of the two atoms to form 
S = si -\- S2. Alkalis have one singlet and one triplet BO 
state, which are described by = and 1, respectively. 
S is then coupled to i to give J = S -\- which is fi- 
nally coupled to the sum of nuclear spins, / = ^1+^2, 
to give the total angular momentum, T. This gives the 
kets \{SI)J^I]TMt)^ which we refer to as the molecular 
basis. The projection Mt is taken along the magnetic 
field axis. 

The entries of K^^^^ depend only on whether a chan- 
nel is of singlet or triplet symmetry, and are given in 
terms of the scattering lengths as,t of the corresponding 
potentials by [22]: 



as,t/a V2 



{BO) 
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• tan(7r/8) 



{BO) 



tan(7r/8) 



(1) 



Here, k[^^^ is the K matrix element of 
the singlet (s) and triplet (t) channels, a = 
2-3/2 [r(3/4)/r(5/4)](2/iC6/n^)^/'^ -g ^^^^ scatter- 
ing length, and r{z) is the gamma function. 

The channel basis by contrast to the molecular ba- 
sis, is constructed in terms of uncoupled atomic states, 
plus the partial wave of their collision. At nonzero mag- 
netic field the projection of the total atomic angular 
momentum /i^2 = ^1,2 + ^1^2 onto the magnetic field 
quantisation axis, mi ^2, is a good quantum number, but 
fi itself is not. The states are then described by the 
kets |aimi, a2m2, ^m^), which we refer to as the Zeeman 
basis. Here, ai^2 = a^b^c... label the energy ordered 
atomic Zeeman states [1]. 

We obtain K^*^ using a frame transformation 2^, 24 1 
from the molecular basis to the Zeeman basis. This 
magnetic- field dependent, unitary transformation can be 
calculated in terms of Clebsch-Gordan coefficients and 
the Wigner 6j and 9j symbols using angular momentum 
algebra [25|. For the case of identical atoms, it is simple 
to incorporate the relevant Bose/Fermi symmetrisation. 
Writing the frame transformation as a matrix U, we have 
K(-) = UK^^^^Ut. 

We next calculate the long range K matrix, K(£^), 
from which the observable scattering properties can be 
extracted. Unlike K^^\ the long range K matrix is en- 
ergy dependent because of the van der Waals tail of the 
potentials. For a multichannel problem with both open 



and closed channels, K(£^) is given by [22[ 
K{E) = -[ZfgiE)-Z,g{E)Kl;^] 



(2) 



where 



K 



is) 
eff 



K'J + K'J [X{E) - K(«,)]-1kW . (3) 



Here, the subscripts 'o' and 'c' refer to the open and 
closed channel blocks of the K matrix. The energy de- 
pendent, diagonal Z matrices are given by the long range 
behavior of the reference functions / and ^, and represent 
the propagation of a wavefunction from small to large 
distances. The diagonal matrix of the bound state phase 
in the closed channels, %(^), is analogously defined [22 1 
and, as shown by the structure of Eq. (|3]), gives rise to 
the resonant enhancement of the collisions. It also allows 
the calculation of bound state energies from the determi- 
nental equation 



det(x(^)-K(«,)) = 0. 



(4) 



The Z and x matrix entries are evaluated for each chan- 
nel I a) at the magnetic field dependent energy E—Ea{B). 
Finally, the S matrix, which we use to extract the ob- 
servable scattering properties, may be calculated from 
S{E) = [l^iK{E)][I-iK{E)]-\ where I is the identity 
matrix. By approximating the short range interactions 
to be energy independent, and using the form V^^{r) 
for the potentials, we have therefore reduced the scatter- 
ing problem to one involving just three parameters: ag, 
at, and Cq- 

Our code is sufficiently fast that it is possible to search 
over the whole (as,at) plane. Given one field at which a 
resonance occurs, and assuming that Cq is known, lines 
through the (ag, at) plane are typically found which result 
in a resonance at this field. Further resonances reduce the 
possible range for each scattering length. For more ac- 
curate determinations, least squares minimization using 
all available resonance locations then gives optimal val- 
ues of as and at- We note that the values thus found can 
be different to those produced by a full coupled channels 
calculation using realistic potentials, and should be inter- 
preted as fit parameters that will in general be close to 
the real scattering lengths. Nonetheless, the scattering 
lengths found can be used for the prediction of unob- 
served resonances. This could make our technique useful 
for investigating new collision partners. 

As an example application of our code we consider 
cold collisions of ^Li and ^^K atoms. The analysis of 
Ref. found = 52.1 ao, at = 63.5 ao, and Cq = 
2322 ^htto- Here, ao = 0.0529177 nm is the Bohr ra- 
dius, and El, = 4.35974 x 10"^^ J. With these input 
parameters, we find s-wave resonances in the aa chan- 
nel at 17.08 mT and 18.13 mT. The locations found differ 
significantly from the measured values of 15.76 mT and 
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TABLE I: Resonance locations Bo and widths AB obtained 
from an MQDT calculation for ^Li-^°K, after varying as and 
at to optimally fit measured resonances. Experimental (Exp) 
and coupled channels (CC) data are from Ref. [10]. Channels 
are labelled by the Zeeman state of each atom, with ai re- 
ferring to ^Li and a2 referring to ^°K, and the total angular 
momentum projection, Mt- 'p' indicates a p wave resonance, 
for which AB is not defined. 







MQDT 


Exp 




CC 


aia2 


Mt 


Bq [mT] AB [mT] 


Bo [mT] 


Bo [mT] AB [mT] 


ha 


-5 


21.36 


0.028 


21.56 


21.56 


0.025 


aa 


-4 


15.93 


0.022 


15.76 


15.82 


0.015 


aa 


-4 


17.01 


0.007 


16.82 


16.82 


0.01 


aa 


-4 


25.80 


P 


24.9 


24.95 


P 


ab 


-3 


0.80 


P 


1.61 


1.05 


P 


ab 


-3 


15.31 


0.042 


14.92 


15.02 


0.028 


ab 


-3 


15.96 


0.016 


15.95 


15.96 


0.045 


ab 


-3 


16.86 


0.005 


16.59 


16.59 


0.0001 


ab 


-3 


26.08 


P 


26.3 


26.20 


P 


ac 


-2 


14.77 


0.044 


14.17 


14.30 


0.036 


ac 


-2 


15.81 


0.056 


15.49 


15.51 


0.081 


ac 


-2 


16.72 


0.005 


16.27 


16.29 


0.060 


ac 


-2 


26.10 


P 


27.1 


27.15 


P 



16.82 mT flO^. However, allowing the singlet and triplet 
scattering lengths to vary, we find optimal agreement for 
tts = 53.17 ao and at = 64.41 ao- These values were ob- 
tained using a least squares minimization, comparing our 
results to the experimentally observed locations of the 
thirteen resonances listed in Table [H Our results, ob- 
tained from only the given Cq and by varying the scat- 
tering lengths, agree well with the more computationally 
intensive coupled channels approach, and with the ex- 
perimental results. The abovementioned resonances, for 
example, are found at 15.93 mT and 17.01 mT, a disagree- 
ment of approximately 1 %. 

The s-wave resonances in the ac channel listed in Ta- 
ble [J are illustrated in Fig. [TJ The two panels show the 
effect of a Feshbach state below and above the threshold 
at which it causes a resonance. The lower panel shows the 
bound state energies as a function of magnetic field. The 
upper panel shows sin^ ^o{E)^ where So{E) is the 5-wave 
scattering phase shift in the ac channel. In this exam- 
ple, there is one open channel and eleven closed channels. 
Consequently, the 5-wave phase shift is linked to the S 
matrix by S{E) = e'^'^^^^^\ Resonances can be recog- 
nised by the sudden change in the value of sin^ So{E). 

As a second example we consider ^^K— ^^Rb scattering. 
In Fig. [2^ we compare an MQDT calculation of 5-wave aa 
channel resonances to three coupled channels calculations 
using two different pairs of singlet and triplet power law 
potentials, and the potentials of Pashov et al. [26[, which 
to the best of our knowledge are the most accurate avail- 
able. The power-law potentials are constructed to have 
the same scattering lengths, Ce, and number of bound 
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FIG. 1: (color online). Bound state energies (lower panel) and 
scattering phase shift sin^ 60 (upper panel) for the ac channel 
of ^Li-^°K, as a function of magnetic field, calculated from 
our three parameter MQDT model. The inset shows a close- 
up of the 15.81 mT resonance, illustrating the narrow region 
within which the Feshbach molecular state has an appreciable 
entrance channel component 



states as the potentials of Ref. [26] . These values for the 
scattering lengths and Cq are also used in the MQDT 
calculation, without varying their values as in the ex- 
ample above, to allow direct comparison. The resulting 
resonance locations illustrate the two main approxima- 
tions of our approach. Firstly, we approximate the short 
range K matrix to be energy independent. The validity 
of this is shown by the close agreement of the MQDT 
result to that of the coupled channels calculation with 
—Cq/v^ + potentials. These quickly converge to 

pure van der Waals potentials as r increases, only hav- 
ing significant C12 contributions at distances where the 
hyperfine energies are small compared to Cq/v^. 

Our second important approximation is the use of a 
pure van der Waals potential. Other dispersion contri- 
butions to the potential, such as the attractive —Cs/r^ 
term, can be important, as shown in Fig. [2^. The 
—Cs/r^ term, taken from Ref. [2§], has a sufficiently long 
range to give deviations from the pure van der Waals re- 
sult. These results are closer to those given by the full 
potential. MQDT using numerical reference functions 
based on a more detailed potential [11] could be more 
accurate, but would be more computationally involved. 

Another way of illustrating the differences between 
the calculations is by studying the near-threshold bound 
states, which we show for the aa channel in Fig. [2)3. The 
difference between the MQDT result and that of the cou- 
pled channels calculation using the potential of Ref. [2^ 
is primarily due to the zero field bound state locations be- 
ing different, which is a consequence of the different sin- 
glet and triplet potentials. The states around —0.4 GHz, 
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neglects the weak magnetic dipole-dipole and second- 
order spin-orbit interactions, which could be included in 
a more general MQDT treatment. 

In conclusion, we have developed a model of Feshbach 
resonances in ultracold alkali gases utilising the separa- 
tion of energy scales suggested by MQDT. A frame trans- 
formation allows us to link short range, where we can ne- 
glect the energy dependence of the interactions, to long 
range, where we can use an approximate potential with 
analytic reference functions. We are therefore able to 
predict and study resonances on the basis of three input 
parameters - ag, at and Cq - in addition to known atomic 
properties. Our approach is significantly less computa- 
tionally challenging than the coupled channels approach 
and, by using the scattering lengths as fit parameters, 
could enable useful predictions to be made while investi- 
gating new collision partners. 

We gratefully acknowledge Bo Gao for supplying us 
with his codes for calculating the Z(£^) and ma- 
trices. P.S.J acknowledges partial support by the U.S. 
Office of Naval Research. 



FIG. 2: (color online). Gomparison of ^°K-^'^Rb aa reso- 
nance locations predicted by different approaches. The up- 
per panel shows scattering length vs magnetic field. The 
MQDT results are given by the approach presented in this 
paper. The other results are given by coupled channels calcu- 
lations with potentials of the form —Cq/t^ + Cx^jr^^ (6-12), 
-Ce/r^ - Cs/r^ + Ci2/r^^ (6-8-12), and the potentials of 
Ref. ^] (full). The full potentials were also used to calcu- 
late the bound state spectrum ('o') which is compared to an 
MQDT calculation ('+') in the lower panel. Bound states are 
numbered according to the vibrational quantum number of 
the singlet and triplet states which give rise to them. Note 
that the two plots have different ranges of magnetic field. 

arising from the most weakly bound {v — —1) level, agree 
to approximately 1%, whereas the more deeply bound 

= — 2 and —3 states deviate more significantly. The 
MQDT bound state energies become less accurate as 
binding energy increases, as non van der Waals contri- 
butions become more important. This implies that reso- 
nances due to more deeply bound levels will be less accu- 
rately reproduced by the MQDT model. The fit proce- 
dure discussed above partially compensates for this. We 
note that the ABM method [lO| has been shown to accu- 
rately predict the locations of ^^K— ^''Rb resonances (27| . 

The MQDT model can give a qualitative indication of 
resonance widths, as shown in Table [D In coupled chan- 
nels calculations, resonance widths are determined pri- 
marily by the difference between the singlet and triplet 
potentials. This is not included in the MQDT model, 
where the only difference between the two is in their scat- 
tering lengths. Lastly, we note that our MQDT approach 
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